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We consider a simple toy model of a regular bouncing universe. The bounce is caused by an 
extra time-like dimension, which leads to a sign flip of the term in the effective four dimensional 
Randall Sundrum-like description. We find a wide class of possible bounces: big bang avoiding ones 
for regular matter content, and big rip avoiding ones for phantom matter. 

Focusing on radiation as the matter content, we discuss the evolution of scalar, vector and tensor 
perturbations. We compute a spectral index of ris = —1 for scalar perturbations and a deep blue 
index for tensor perturbations after invoking vacuum initial conditions, ruling out such a model 
as a realistic one. We also find that the spectrum (evaluated at Hubble crossing) is sensitive to 
the bounce. We conclude that it is challenging, but not impossible, for cyclic/ekpyrotic models to 
succeed, if one can And a regularized version. 

PACS numbers: 04.50.-|-h,98.80.-k,11.25.Wx,98.80.Es 


I. INTRODUCTION 

Recently, bouncing models of the universe were re¬ 
vived in the framework of string cosmology Q, since 
it is possible to generate a scale invariant spectrum 
of density fluctuations in the pre-bounce phase. Spe¬ 
cific realizations are e.g. cyclic/ekpyrotic models of the 
universe i II i i i or the pre-big bang scenario 
The main problem for such models to suc¬ 
ceed lies in the fact that the bounce itself is often sin¬ 
gular in the models studied so f ar llfl. T his makes 
matching conditions a necessity |llL ll2L IlIlL Q which 
are, unfortunately, rather ambiguous. In addition, fluc¬ 
tuations often become non perturbative near a singu¬ 
larity M In the few toy models of regular bouncing 
cosmologies known in the literature, the actual bounce 
has a stron g impac t on the evolution of perturbations 

El El Eiiii M El m m m m [23. what is 

more, the models are often technically challenging at the 
perturbative level and unfortunately inconsistent meth¬ 
ods were proposed in the literature |^. Consequently, 
bouncing models have not been taken seriously as alter¬ 
natives to the inflationary paradigm 

In this article we will construct a simple regular toy 
model, for which one can compute analytically how per¬ 
turbations evolve through the bounce. We find a red 
spectral index for scalar perturbations that is ruled out 
by observations |^. Furthermore, the detailed mecha¬ 
nism providing the bounce has a strong impact on the 
spectrum. 

In order to generate a bounce within pure four dimen¬ 
sional general relativity, one has to violate various energy 
conditions. The usual mechanism in the literature in¬ 
volves the addition of an extra exotic matter field that be- 
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comes dominant during the bounce 

(see also for an a’ regularized bounce. |3lL 1 ^ for 
a non-singular cosmology achieved by high er derivative 
modification of general relativity, or for a bounce 

induced by a general-covariant, T-duality-invariant, non¬ 
local dilaton potential). The presence of a second matter 
held, and perturbations in it, is one physical origin for 
the sensitivity to the type of bounce one is considering. 

We will not add an additional matter held, but use an 
ingredient motivated by string theory: extra dimensions. 
If we assume one additional dimension, the effective four 
dimensional Friedmann equation will show corrections 
proportional to at high energy densities [H3- If the 
extra dimension is space-like, we have the usual Randall 
Sundrum setup . In that case, p and p^ have the 

same sign, and no bounce will occur. 

However, if we assume the additional dimension to be 
time-like, this sign will flip and thus it will cause a bounce 
when p^ becomes dominant . We are aware 

of complexities associated with an extra time-like dimen¬ 
sion; for instance, existence of tachyonic modes and pos¬ 
sible violation of causality, or pos sible appearance of neg¬ 
ative norm states (see Ell for more discussion 

regarding these issues). We will not touch on those issues, 
but take the modification of the Friedmann equations as 
a simple toy model that could also arise by other means. 
The outline of the article is as follows: first, we work 
out the background solution based on the above idea in 
section im We find a whole class of simple, analytically 
known, bouncing cosmologies, either big bang avoiding 
ones, for usual matter content, or big rip avoiding ones, 
for phantom matter. We will then focus on one specific 
big bang avoiding bounce by considering radiation as the 
matter content - this seems to be the most conservative 
choice to us. 

In Section imi we discuss scalar perturbations and com¬ 
pute analytically the spectrum of fluctuations in the 
post bounce era, for Bunch-Davis vacuum initial con¬ 
dition in the pre-bounce era. We find a red spectrum 
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with a spectral index of ns = —1, thus ruling out this 
model as a realistic one. We also find that the back¬ 
ground cosmology around the bounce dictates the de¬ 
formation of the spectrum through the transition. In 
Section liin^ we conclude that it is challenging for the 
bounce in Cyclic/Ekpyrotic models of the universe to 
leave the spectrum unaffected. The only way to check 
if this is the case is by finding a regularized version of 
the proposed scenarios and following the perturbations 
explicitly through the bounce. 

In Section im we discuss vector perturbations and find 
that they remain perturbative during the bounce. Fi¬ 
nally in Section m we derive analytically the evolution 
of gravitational waves as they pass through the bounce. 
We find a blue spectrum of gravitational waves in the 
post bounce era on large scales, and an amplitude that 
depends on the details of the bounce. 

II. BACKGROUND 


We may rephrase the quadratic corrections in terms of a 
second ideal fluid with 


P- ■■= ^ (2p + p) , (7) 

where A = The ’’two fluids” are of course related 

and represent only one degree of freedom. If we have an 
equation of state between := p and p+ := p, that 
is p+ = wp+, we get p_ = W-p- with w- = 2w + 1. 
It will become useful to write w = (n — 3)/3 so that 
w- = 2n/3 — 1. As a consequence of the equations of 
state we get 

p+{t) = ra{t)-'^, (8) 

p_{t) = . (9) 

The Friedmann equation reads now 


We use a metric with negative signature, scale factor 
a{t), cosmic time t and, for simplicity, we work with a 
flat universe, so that 




and has the solution 


( 10 ) 


= dt^ — a^Sijdx^dx^ . ( 1 ) 


Considering a modified version of the Randall Sundrum 
(RS) scenario [^1^ . or iginating e.g. by having an extra 
time-like dimension | 3 !tL I3B| we have the modified Fin- 
stein equations 


G — K -j- K Sfj^u , 


( 2 ) 


where = Sir/Mp, = Sir/Mp and Mp is the fun¬ 
damental 5-dimensional Planck mass. For simplicity and 
since we are primarily interested in the high energy mod¬ 
ifications due to Sp,^, we fine tuned the model such that 
the effective four dimensional cosmological constant van¬ 
ishes and we neglected the projected Weyl tensor. 

Spi, is quadratic in the energy momentum tensor Tp^, 
and is given by ^2) 


Sa. = 


^ rjia rp _ 

T ii.ij . -L n.r\-L , 
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(3) 


We have a different sign in front of Spi, in m compared 
to the usual RS setup [33. This modification will yield 
a class of non singular bounces we shall examine below. 
If we consider an ideal fluid with 

m p".), (4) 

the quadratic term becomes 





0 

^p{p + 2p)5] ) 


(5) 


a{t) 


r 





( 11 ) 


for n 7 ^ 0. Thus for any positive n we get a regular 
bounce at t = 0 with minimal scale factor 


ao = (r/2A)i/" , (12) 

e.g. for dust (n = 3), radiation {n = 4), stiff (holo¬ 
graphic) matter (n = 6) etc. Negative n corresponds 
to matter that violates the null energy condition (NFC). 
However, there is no big rip, but a bounce at the maxi¬ 
mal scale factor oq. Such a phantom bounce was recently 
proposed in [^, but no analytic solution was given. If we 
also introduce the characteristic time scale of the bounce 



we may write 

a{t) = ao 1^1 -f 


(13) 


(14) 


The relation between cosmic time t and conformal time 
r] is given by 



(15) 

(16) 


where F^{x) is a generalized hypergeometric function, 
but we will not make use of this relation in the following 
analysis. 
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Let us turn our attention to perturbations and exam¬ 
ine how they evolve through a big bang avoiding bounce 
with minimal scale factor. This is the most interesting 
case 1^ , because it could be seen as a simple toy model 
for a regularized version of the bounce occurring in the 
cyclic scenario. It is a toy model, since we remain at the 
four dimensional effective description throughout. We 
will comment on this issue in section Fill PI 

In the following, we have to be careful, since the so 
called “sound speed” cf = (p+ —p_)/(/?+ — /9_) diverges 
near the bounce To be specific, for our background 
solutions one gets 

2 (a/aor(n/3-l)-(4n/3-2) 

'• = - 


and from perturbing m we get 

.2 / 




Po 


{l + w)V, 


(23) 


6 \-{l + w)Vi [(1-f 2w)^-h s](5* 

With Ti. = a'/a the perturbed Einstein equations read 

(25) 


- 3H(H$ -k $') = —K- 


-b + {2-H' + = y kVo 


-b s - [(1 -b 2w)^ -b s] 


m + ^) (i - . (26) 


so that Cg diverges at a” = 2^/”aQ, that is at t = ±to- 
As a consequence, the usual equations governing the evo¬ 
lution of perturbations that involve can not be used 
near the bounce. Therefore, we will derive the relevant 
equations of motion from first principles in the next sec¬ 
tion. 


III. SCALAR PERTURBATIONS 


For the time being we will work with conformal time rj. 
The most general perturbed metric involving only scalar 
metric perturbations in the longitudinal gauge is given 

by lilllllil 

ds^ = [{1 + 2^)dri^ — {1 — 2'^)Sijdx’'dx^] -(18) 

Note that the two Bardeen potentials agree with the 
gauge invariant scalar metric perturbations - there is no 
residual gauge freedom. 

The perturbed energy momentum tensor is given by 

Ei 


((5^) 


f -Sp {po+Po)Vi\ 

V -{po+Po)Vi - cr-,ij j ’ 


(19) 


where po = P+ and po = wpo denote background quanti¬ 
ties given by and V is the velocity potential so that 


Note that the right hand side of (1^ will become zero 
at dP = r/A. Henceforth, one might expect difficulties 
in deriving a single second order differential equation for 
$. A method to deal with this problem, in the case of 
adiabatic perturbation or s = O^was advocated in E3 
and subsequently used e.g. in l3C]| and the previ¬ 
ous versions of the article: one splits <I> in two compo¬ 
nents, each of which satisfies a regular second order dif¬ 
ferential equation. However, this method is inconsistent 
with other constraint equations, like the fluid conserva¬ 
tion equations in the case of true two fluid models, as was 
explicitly shown in , or in our model the fact that the 
density of the second fluid is related to the first fluid. 
There, an alternative method was introduced, which we 
shall employ in the following. 

Since the universe was radiation dominated before the 
period of recombination, let us restrict ourselves to n = 4 
such that 

a{x) = ao [l + , (27) 

Where, a; := -A is the rescaled time. Now, let us combine 
and to 

0 = (l-g2u;)^) _ 37^(7^$ + $')] 

- (i - [$" + + {2n' -b n^)^] ,(28) 


= ( W ) ■ (20) 

Considering no anisotropic stress, that is a-^ij = 0, the 
off-diagonal Einstein equations yield 


$ = 4/, 


( 21 ) 


where we used po = since n = 4 and w = 1/3 for 

radiation. Converting this equation to rescaled cosmic 
time X, making use of the background solution 12711 for 
a (a:) and introducing the rescaled momentum 


ao 


(29) 


so that only one scalar metric degree of freedom is left. 
Considering 6p = w5p + tSS, where SS is the entropy 
perturbation, and defining ^ := Sp/po as well as s := 
tSS/po we have 


(sn) = PO 


-e a+w)v,\ 

-{l + w)Vi (wC + s)<5) J 


( 22 ) 


yields 


1 - 


l-ba;2 
(1 -b a;2)2 




1 X 


21 +x^ 
1 10 1 


5- 


18 


3 3 l-ba;2 


l-ba;2 


vT" 






(30) 

0 , 
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where is a Fourier-mode of $ and dot denotes a deriva¬ 
tive with respect to x. 

Before we go ahead and solve this equation in differ¬ 
ent regimes, let us have a look at the problematic region 
around a: = ±1, the boundaries of the region where the 
Null Energy Condition is violated. The nontrivial ques¬ 
tion is whether <i>fe remains continuous and smooth near 
these points or not. This question has been the object of 
study in |25| , which provides analytic solutions and a de¬ 
tailed discussion in the case of adiabatic perturbations. 
Consequently, we will expand near these points to 
find the two independent analytic solutions, as proposed 
in m. 


A. Approximate analytic solution 

Our final goal is to compute the spectrum 

7?^, := |$2| (31) 

long after the bounce {x 3> 1), while specifying the ini¬ 
tial conditions for each Fourier mode well before the 
bounce (x <C — 1). To achieve that, we will find ap¬ 
proximate analytic solutions within different regions and 
match them together. In this way, we will be able to 
compute analytically the spectral index Ug. 

Let us start in the limit x —>■ —oo. In that case one 
can introduce the gauge invariant variable Vk in terms of 
which the action takes the simple form of a scalar field 
(see for a detailed discussion) so that the equation 
of motion reads 


+ 

{clk'^-—\vk = 0, 

(32) 

\ z J 

with Cg = 1/3 and 

jl 

(33) 


P := T-i^-n'. 

(34) 


in that regime. If we use the background solution in the 
limit a; <C — 1 we get 




33/4 j 

2 k^/'^x 


V3 1 \ 

2 j 


exp 



(39) 


as an approximate analytical solution for $. The overall 
scale is given by 


a := 


Ip^/t^ 
3/2 ■ 


(40) 


which is a free parameter in our model, subject to mild 
constraints, e.g. to should be larger than the Planck time, 
but smaller than the time of nucleosynthesis. However, 
we will not need to specify its value if we scale the spec¬ 
trum appropriately. 

Similarly, for a: 3> 1 the general solution of 


5 • 

‘bfe H-H-= 0 

9t 3t 


can be written as 

^ i-^ky/x /. 1 

= —e (j- ) 

X 2 ky/x 

X 2 kyCv 

where the Ci are constants. For 

1 < bl < 




(41) 


(42) 


(43) 


we can neglect the fc-dependence of pnji . that is we can 
use 

= (44) 


f 2 / 

1 •• 1 , 

( 18\ 

1 - ^ 

+ — 

5- ^ 

V ) 

' 2x ' 

1 x^ J 


Since the universe is radiation dominated, z" vanishes 
and we can solve (1221 




(35) 


where we imposed quantum mechanical initial conditions 
for Wfc, that is 


Vki.'nQ) = k (36) 

v'kiVo) = ik^l'^N (37) 


and the normalization condition NM* + N*M = 2. 
The variable Vk is related to 4)^ via 


which is solved by 

= C3(x2-2)2a;-(7-'/4T)/4^ni^^2/2) 

Ci{x^ - 2)2a:-(^+'^)/^F”;(xV2), 

where 

ni = 

di = 

742 = 


9-hV4T is-hViT 


1-f 




15- 9- ViT 


(45) 

(46) 

(47) 

(48) 





/3 



(38) 


C?2 = 


~T' 


(49) 
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and FJl{z) is the generalized hypergeometric function. 

Close to a; = ±1 we can Taylor-expand so that we 
get 


e$fc + 





1 

4 



(50) 


where we introduced 


X =: ±1 + e . (51) 

We will now follow to solve this equation near x = 
± 1 : if remains bounded, we get one independent so¬ 
lution from which we can deduce the other one by means 
of the Wronskian method. In this way, we arrive at 

chfc = Cs (l - e) + Cee'. (52) 


Note that it would be justified to neglect the k depen¬ 
dence in this solution. For the detailed reasoning as to 
why this solution is appropriate we refer the reader to 

Ill- 

Last but not least, we have to find a solution for the 
actual bounce, that is for |x| <C 1 and fc <C 1. The 
relevant equation of motion in this regime is 


-f yX$fc -f =0, 


(53) 


which is solved by 

-13a:^/8 

^ (C7lT™(13xV7) + C8M™(13 xV7)) (54) 

V X 

with m = —9/52, n = 1/4 and W,M Whittakers func¬ 
tions. Note that the solution is well behaved throughout 
the actual bounce. 

All that is left to do is to match all the solutions 
smoothly. To be specific one has to use for —oo < 
X < —X, for —X < X < —1 — e, for —1 — e < 
X < —1 + e, it^ for —1 -I- e < x < 1 — e, (1^ for 
1 — e<x<l-|-e, lHIHl for 1 -I- e < X < X and finally 
for X < X < 1/(4A;) with some small e. This straightfor¬ 
ward but tedious calculation results in a spectral index 
of 


Us = -1, (55) 

which is independent of the choice of e This in¬ 

dex is of course in conflict with the observed nearly scale 
invariant spectrum. 

A few words regarding the validity of the matching pro¬ 
cedure might be in order. Lets first consider the match¬ 
ing at X = ±x: To either consider the term oc or the 
term stemming from H -|- 2H^ in (EOl) and match at the 
point where both are equal is indeed a standard proce¬ 
dure that was employed and tested (e.g. numerically) at 
various instances. For example, in the framework of a 


bouncing universe it was used in M and verified numer¬ 
ically in the limit fc <C 1, [b^. One can simply check the 
validity of this matching for our case since the solutions 
of (HHll as well as equations ll^ and ( 1 ^ . around |x| ~ x 
behave as 

(56) 

where A^’s and B^’s are constant. This ensures a smooth 
matching of these solutions at these points. Also, one 
could successfully use the same matching procedure in 
the study of fluctuations in an inflationary universes, 
even though a more elegant procedure is available (see 
e.g. Chapter 13 of ^ detailed discussion). The 

remaining matchings occur close to ±1 and do not influ¬ 
ence the spectral index as long as we are confident that 
the solutions are regular (as is the case for our scenario), 
since all equations dictating the matching conditions are 
independent of k in the limit of fc ^ 1 . 

As a consequence, we will be able to give a simplified 
analytical argument in the next section yielding the same 
spectral index. Thereafter, we com par e our result with 
related studies in the literature like . 


B. Simplified matching procedure 

We would like to provide some simple insight into our 
result of Ug = — 1 from the previous section. As we men¬ 
tioned earlier, it can be shown that the solutions of the 
corresponding equations around |x| ^ x have the follow¬ 
ing asymptotic x dependence 


Bi 

^k-Ai + ^ 3^2 

X —x~ , 

(57) 

= As + 

X ^ — x^ , 

(58) 

B-i 

— A^+ ^3^3 

X ^ +x~ , 

(59) 

. i?4 

- ^4 + ^3/2 

X ^ +x~^ . 

(60) 

This makes an analytic tracking of the calculation very 
easy: on the one hand, the smooth matching conditions 

are satisfied by requiring 



Ai = As, 

Bi = B 2 

(61) 

As = A 4 , 

B 3 = B 4 , 

(62) 


and on the other hand, the proper k dependence of Ai 
and Bi can be enforced through the the appropriate ini¬ 
tial conditions leading to ra . Taylor expanding in 
the limit <C 1 results in 

^ ^2 ^3^^ “ “ A:(-x)3/2 ^ 
so that we have 

Ai oc (64) 

Bi oc . (65) 
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Next, note that the transfer functions relating A 2 and 
B 2 to the post-bounce coefficients A 3 and B 3 have to 
be independent of fc, since the equations governing the 
regime in between —x and x are independent of k. In 
fact, our equation have regular solutions throughout this 
region so that it is reasonable to deduce 

^3 = £1^2 + £2^2 1 ( 66 ) 

B 3 = 1 - 1 A 2 L 2 B 2 , (67) 

where ei, £ 2 , (-i and i 2 have to be some constant numbers. 

Consequently the total transfer functions will be 

A4 = £1^1 + £27?i , (68) 

B 4 = ^.1^1 tBi. (69) 

Finally, making use of the above equations and A 4 ^ Bi 
for small k according to lIHlll . we can calculate the power 
spectrum for A 4 , the non-decaying mode of (HU 

Vk ■= k^ \A^\ (X k^ \Bf\ cx k~^ . (70) 

Henceforth, we conclude a spectral index of 

ns = -l. (71) 


C. Numerical treatment 

In order to provide a quick check of our analytic ar¬ 
guments, we integrated equation (TTnil with a 4/5th-order 
Runge-Kutta method, as implemented in MAPLE v9. 
Starting the numerical treatment at —x with initial con¬ 
ditions given by and evaluating the final spectrum 
at X = l/(4fc) yields a spectral index of ris = —1, Fig^ 
The computations in MAPLE were done with 120 dig¬ 
its accuracy; nevertheless, it should be noted that the 
amplitude can not be recovered properly by the simple 
numerical method used, due to the non-trivial nature of 
the equation of motion near the bounce. However, this 
does not effect the spectral index since the actual bounce 
is /c-independent, as can be seen in 1S3). 


D. Interpretation 

Our result of the previous section clearly shows that for 
the specific bouncing cosmology considered, both modes 
(the initial constant one and the growing one) will pass 
through the bounce, but it is the mode with the redder 
spectrum, i.e. the growing mode for the Bunch-Davis 
vacuum, that dictates the final shape of the post-bounce 
spectrum. Henceforth, we observe mode mixing in a 
model with only one degree of freedom and no spatial 
curvature. One has to be more careful to draw such con¬ 
clusion for more general cases involving two different flu¬ 
ids (constrained by an adiabatic condition in other back¬ 
grounds), since the matching procedure at ±i is not that 
simple in general. To be more specific, the feature of the 



FIG. 1: The numerical integration of yields the spectrum 
iinii evaluated at a; = l/(4fc) (dots); the solid line corresponds 
to a spectral index of ris = —1, in agreement with the analytic 
result 11551 . Higher order effects in k are accounted for by 
letting k cover a broad range from 1 down to 10“^°. 


bounce which is crucial for determining the approximate 
behavior of the out coming spectrum, besides knowing 
the existence and regularity of the solutions throughout 
the bounce region, is the correction to the equation of 
motion (Ell due to the second fluid. This is the region 
where modes just turned non oscillatory and corrections 
due to the bounce have to be considered, because they 
are crucial for the solution. 

Recently, a wide class of bouncing cosmologies with 
two independent fluids was discussed in by Bozza 

and Veneziano [^ . One might think that, in their no¬ 
tation, our class of bouncing universes corresponds to 
q; = 0, = (n — 3)/3, cl = 2n/3 — 1 = 2 c^ -f 1 and 
n = 4. This would put it into region C of and a 
spectral index of 


Hs 


4(n-3) 
n — 1 


3 


(72) 

(73) 


would result according to [^. However the model of 
differs from our bouncing universe, since a = 0 only 
enforces the adiabatic condition on each fluid separately 
and while no intrinsic entropy modes for each fluid are 
considered, two independent initial conditions for each of 
these fluids are assumed. In our model the ’’two fluids” 
are indeed related, so that the total pressure becomes a 
function of the total energy (as a consequence only one 
initial condition can be specified). In other words, we 
have enforced the adiabatic condition in its strong inter¬ 
pretation. It should also be noted, the approximation 
scheme employed there, differs considerably from ours, 
involving an expansion in k and a dimensional argument 
for the contribution of the bounce (no explicit solution for 
the bounce was used). Henceforth, the results of [23, [^ 
can not be applied to our model. 

A few words rerarding Ekpyrotic/Cyclic scenarios are 
needed: in y, 13 0 there is a mechanism present that 
generates a scale invariant spectrum before the bounce 
13 . It was shown in that this spectrum will not 
get transferred to the post-bounce era, if perturbations 
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evolve independently of the details of the bouncing phase. 
Furthermore, the main conclusion of was, that a 

smooth bounce in four dimensions is not able to gen¬ 
erate a scale invariant spectrum via the mode mixing 
technique. However, all these arguments were specific to 
bouncing cosmologies in four dimensional gravity. So far, 
the effective four dimensional description breaks down 
for all proposed cyclic models [g near the bounce, and 
the bounce itself is singular i M- Thus, it is of 
prime interest to search for a fully regularized bounce in 
cyclic/ekpyrotic models of the universe, so that one can 
compute the spectrum of perturbations explicitly in 
a full five dimensional setting - see e.g. for a recent 
proposal and for argument that mode mixing can 
occur in the full five dimensional setup. Our results con¬ 
firm that it is challenging for such a model to succeed, 
but not impossible. 


IV. VECTOR PERTURBATIONS 

Vector perturbations (VP) have recently caught atten¬ 
tion in the context of bouncing cosmologies, due to their 
growing nature in a contracting universe [bill l54| . In the 
context of a regularized bounce that we are discussing, 
VP remain finite and well behaved, as we shall see now. 

The most general perturbed metric including only VP 
is given by 0, ^ 


(%.) = -a" ( , (74) 

where the vectors S and F are divergenceless, that is 
S\i = 0 and FV = 0. 

A gauge invariant VP can be defined as 

CT* = S^ + F^'. (75) 


Newtonian gauge. For simplicity we assume the absence 
of anisotropic stress, tt* = 0 , so that 


iSS^.) = qPo 


0 ip+p)V^ 

-{p + p){v^ + S^) 0 


(79) 


The resulting equations of motion read 
1 


2 a 2 


= -f 


p„(l + .B)(l-y)y, (80) 

“ = (“’(*!■ +A)) ■ ( 81 ) 


Equation m is easily intergrated for each Fourier mode, 
yielding 




(82) 


where (7)^^ is time independent. Since a > 0, there is no 
divergent mode present. Equation p1|l is an algebraic 
equation for each Fourier mode of the velocity perturba¬ 
tion U*. Since (1 — po/A) becomes zero shortly before 
and after the bounce, it seems like U® has to diverge be¬ 
fore the bounce - however, all that this is telling us is 
that V* is not the right variable to focus on: we should 
focus on the combination of U® and po that appears in 
the perturbation of the effective total energy momentum 
tensor, that is po(l + ?i') (1 — Po/A) Eb This combination 
is well behaved, and in fact one can check that it is al¬ 
ways smaller than either p+ — p_ or p+ — p_ if we impose 
appropriately small initial conditions for . Thus the 
perturbative treatment of VP is consistent. 


V. TENSOR PERTURBATIONS 


The most general perturbation of the energy momen¬ 
tum tensor including only VP is given by 

- f ® (po+Po)E* \ 

^ - [-(p,+p,)(V^ + S^) _po(7rV+7r^;j J ’ 

(76) 

where tt® and E* are divergenceless. Furthermore the 
perturbation in the 4-velocity is related to E* via 

{Sun = (fo) ■ (77) 

Gauge invariant quantities are given by 

0 * = V^-F^' (78) 

and 7rb 

From now on we work in Newtonian gauge, where 
F* = 0 so that cr coincides with S and 9 with E. Note 
that there is no residual gauge freedom after going to 


In this section we calculate the amplitude and spec¬ 
tral index of gravity waves |68|. corresponding to tensor 
perturbations of the metric, in a radiation dominated 
universe. Because the complete five dimensional calcula¬ 
tion is very complicated we will continue using the the¬ 
ory of linear perturbations in a four dimensional effective 
field theory. We include only the corrections oc due 
to fifth dimension, entering through the evolution of the 
background. 

In the linear theory of cosmological perturbations, ten¬ 
sor pert urbations can be added to the background metric 
via |46l| 

= n[dn — {Sij + hij)dx^dx^ , (83) 

where hij has to satisfy 

n = h}d=0. (84) 

Thus hij has only two degrees of freedom which corre¬ 
spond to the two polarizations of gravity waves. One can 




FIG. 2: V(a;) sets a scale for the behavior of each mode, it 

rises up to | around the bounce but falls off very quickly as 
X ±oo. 


show (see that the amplitude h for each of these 
polarizations has to obey 

M = 0, (85) 

a 

where /x = ah. Note that gravity waves do not couple 
to energy or pressure and that they are gauge invariant 
from the start. 

If the second term dominates over fj,a"/a in the above 
equation /i will exhibit oscillatory solutions, and if the 
gradient term is negligible we will have /i ~ a, or, in 
other words, h is almost constant. Using the standard 
Fourier decomposition we get for each Fourier mode 

hk + Mfe = 0. (86) 


Therefore, a”/a sets a scale for the behavior of each 
mode. This is very similar to the role of the Hubble 
radius for inflationary scenarios, but rather than being 
almost constant, a”/a is more like a potential barrier in 
quantum mechanics: 


(^ ^ 
\ao/ a 2 (1 + a:^)^/^ 


(87) 


As Fig. 0 shows, V{x) rises up to ^ around the bounce 
but falls off very quickly as x —> ±oo. 

Thus, for i, the solutions of equation are 

oscillatory. This implies that we can set our initial condi¬ 
tions according to the Bunch-Davis vacuum at the initial 
time r]i. We get 

= ( 88 ) 


spectrum of perturbations and the spectral index for 
these modes to be 


Vk ■■= 


nt = 


2k^ I 2 

2 . 


2 e 


oc A:"* , 


(89) 

(90) 


implying that, for ^ we are simply left with 
vacuum fluctuations with an ever decaying amplitude. 

Let us now turn our attention to the more interesting 
case: k"^ 5 . As Fig. 0 demonstrates for these modes, 

we have transitions from the regime k^ ^ ^( 2 :) (i.e. k^ ^ 
a"/a) to the regime k"^ V(x) at the transition points 
±a:j^(±? 7 j^). Here x^^ iVtr) i® taken to be positive and 
using m we get 


X,, ~ 2-1/3 fc-2/3 . (91) 


For |x| Xj^ (1 7 ?I ? 7 j^) we can simply approximate the 

solutions of equation p7|l to be oscillatory 

(92) 

Mv) = ^ , V » V.. (93) 

where Ai and Bi are constants set by initial conditions. 
In the case of the Bunch-Davis vacuum we can set Ai = 1 
and Hi = 0. We will derive A 3 and B 3 by matching the 
solutions smoothly at the transition points. 

But first, we need to estimate the solutions for k"^ 
V(x) (i.e. k'^ <C a"/a) that is for |x| <C x^^. Unfortu¬ 
nately, unlike the case of most potential barriers in Quan¬ 
tum Mechanics, we can not use the WKB approximation 
to calculate the outgoing spectrum, since the condition 


djk^ - a"/a)/dr, 

{k^-a"/ay 


(94) 


is not satisfied in this regime. Fortunately we can still 
approximate the solution of equation (IHSll perturbatevely 
in orders of 


hkiv) 

0 ( 77 ) 


A 2 


l-fc2 




-I-H 2 



k^ 



(95) 



-f • • • . 


For 1 77 1 ^ rj^^, since 


Jo a? 


2 .m% 


ao\x 


1/2 


aok 1/3 


(96) 

(97) 





(98) 


long before the bounce and it will remain basically un¬ 
changed throughout the bounce. We easily compute the 


and also 


0 


0 
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one can further approximate the result of equation lEH) 
to 


/ife(x) = A 2 ao\x\^^‘^ + 52 % —[2.6|x|^/^ - 2]. (99) 


On X 


Notice that = 0 for the above solution, since 

^200 


dr] 


±- 


to 


B 2 

2.60— 

ao 


Const. 


( 100 ) 


at I 77 I « 77 j^. Therefore, /ifc can be smoothly matched to 
our solutions from equations © and at the transi¬ 
tion points iXj,, by equating ]ik and ]i'^. We obtain the 
transfer function by relating the coefficients Ai and Bi 
to A 3 and B 3 . These satisfy the relations 

^3 ^ ifc-i-h2.32fc-i/3]^i (101) 


-f6-*'=’'-[0.65 ik-^]Bi, 


B 3 ~ 


Jkr)i 


[-0.65 ik-^]Ai (102) 

j^(3-ik(2ritr-+^i) [_o.65 ik-^ + 2.32k-^^^]Bi. 


The above result is valid for general initial conditions. We 
first note that, to lowest order in k~^, the constant mode 
in the pre-collapse phase (growing mode for h), matches 
onto the growing mode in the post-collapse phase (con¬ 
stant mode for h), without any change in the spectral 
shape. As we mentioned before, for a Bunch-Davis vac¬ 
uum one simply takes Ai = 1 and Bi = 0. The resulting 
power-spectrum for such a case is 


Vk = 13.44 




„-i 


sin^ {k{r] - ??^J) . (103) 


The above spectrum is oscillatory and, for k{r]—ri^^) ^ 1, 
it is oscillating so fast that one could conclude an almost 
flat {rit = 0 ) effective spectrum, with a decaying am¬ 
plitude. However, if k(r] — ry^^) <?; 1 or in other words 
1 X <C 1/(4A;^), which is the case for super-Hubble 
wavelengths, the above equation simplifies to 


Thus, we get a non-decaying mode which has a blue 
spectrum (rit = 2). Notice, the amplitude of the power 
spectrum was dictated directly by the scales and details 
of the bounce. This result is in agreement with the con¬ 
clusions of 0 , where a scale factor somewhat similar to 
our background scale factor was used. 


the evolution of perturbations during singular bouncing 
cosmologies in string cosmology, occurring e.g. in the 
cyclic scenario. The bounce is caused by the presence 
of an extra time-like dimension, which introduces correc¬ 
tions to the effective four dimensional Einstein equations 
at large energy densities. If matter respecting the null 
energy condition is present, a big bang avoiding bounce 
results, and for phantom matter a big rip avoiding bounce 
emerges. We find analytic expressions for the scale factor 
in all cases. 

Having the Cyclic scenario in mind, we focused on a 
specific big bang avoiding bounce with radiation as the 
only matter component. We then discussed linear scalar, 
vector and tensor perturbations. 

For scalar perturbations we compute analytically a 
spectral index of ris = — 1 , thus ruling out such a model 
as a realistic candidate. We also find the spectrum to be 
sensitive to the bounce region, in agreement to the ma¬ 
jority of regular toy models discussed in the literature. 
We conclude that it is challenging, but not impossible, 
for Cyclic/Ekpyrotic models to succeed, if one can find a 
regularized version. 

Next, we discussed vector perturbations which are 
known to be problematic in bouncing cosmologies, due 
to their growing nature. We checked explicitly that they 
remain perturbative and small if compared to the back¬ 
ground energy density and pressure. 

We concluded with a discussion of gravitational waves. 
The transfer function was computed analytically and we 
showed that the shape of the spectrum is unchanged in 
general. In the special case of vacuum initial conditions, 
the spectral index in post bounce era for super-Hubble 
wavelengths is blue. However, the amplitude turned out 
to be sensitive to the details of our model. 

To summarize, having an explicit realization of e.g. 
the Cyclic model of the universe that features a regular 
bounce it might be possible to produce a scale invariant 
spectrum. Henceforth, the challenge for the Cyclic model 
is twofold: 

1. To find a convincing mechanism that regularizes 
the bounce. 

2. To deal with the perturbations in a full five dimen¬ 
sional setup and to take into account the presence 
of an inhomogeneous bulk as inevitable by the pres¬ 
ence of branes. 

Preliminary results show that, by incor por ating ideas of 
string/brane gas cosmology (see for recent re¬ 

views), a regularized bounce is possible (in preparation). 


VI. CONCLUSION 

We constructed a wide class of regular bouncing uni¬ 
verses, with a smooth transition from contraction to ex¬ 
pansion. This was motivated by open questions regarding 
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